Solving the cosmic lithium problems with primordial late-decaying particles 
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We investigate the modifications to predictions for the abundances of light elements from standard 
Big-Bang nucleosynthesis when exotic late-decaying particles with lifetimes exceeding ~ 1 sec are 
prominent in the early Universe. Utilising a model-independent analysis of the properties of these 
long-lived particles, we identify the parameter space associated with models that are consistent with 
all observational data and hence resolve the much discussed discrepancies between observations and 
theoretical predictions for the abundances of 7 Li and 6 Li. 

PACS numbers: 26.35.+C, 98.80.Cq, 98.80.Ft 



I. INTRODUCTION 

Standard Big-Bang nucleosynthesis (SBBN) is one of 
the most reliable and farthest reaching probes of early 
Universe cosmology, being based on the rigorously tested 
Standard Model of particle physics as well as basic prin- 
ciples of nuclear physics. Augmenting these principles 
with experimental data concerning nuclear reactions, we 
can precisely estimate the relative abundances of light 
elements (ALEs), particularly D, 3 He, 4 He, 6 Li, and 7 Li 
(relative to H) at the end of the "first three minutes" after 
the Big-Bang, as a function of the baryon-to-photon ra- 
tio r). Consequently, SBBN has been instrumental in the 
realisation that baryonic matter constitutes only a small 
proportion of the total energy density of the Universe, 
hence providing further supporting evidence for the exis- 
tence of dark matter. Utilising the value of rj measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP), 
rj WMAP = (6.10 ±0.21) x 10~ 10 at the 68% C.L. Q, the 
majority of theoretical predictions from SBBN are in ex- 
cellent agreement with observational data. This is truly 
remarkable considering that the ALEs span many orders 
of magnitude from 4 He/H- 0.08 down to 7 Li/H- 10 -u . 

Despite the great success of SBBN, it has been noted 
that the prediction for the ratio of 7 Li to H, —9.476 < 
log 10 ( 7 Li/H) < -9.322 (obtained using rj WMAF ), does 
not agree with current observations. Recently Bonifacio 
et al. 0] reported that logio( 7 Li/H) obs = -9.90 ± 0.09, 
confirming the results of Ryan et al. [|[ logio( 7 Li/H) b s = 
—9.91 ± 0.10. Even the less stringent observational re- 
sult -9.63 ± 0.06 of Melendez and Ramirez (MR) ex- 
ceeds the SBBN prediction at 2 a. This is the " 7 Li prob- 
lem". It has been often argued that the 7 Li abundance 



would be smaller than the SBBN value due to depletion 
in stars [H, [||. However, a quite uniform 7 Li abundance 
(the Spite plateau) observed in metal-poor stars is some- 
what difficult to explain with such stellar depletion and, 
moreover, the recent detections by [7| of more fragile iso- 
tope 6 Li in some of these stars provide evidence against 
stellar depletion. 

The observations of [3] that some metal-poor stars have 
the isotopic ratio 6 Li/ 7 Li of a few percent not only sharp- 
ened the 7 Li problem but also elucidated the " 6 Li prob- 
lem". Since the 6 Li abundance predicted in SBBN is 
very small ( 6 Li/ 7 Li- 3.3 x 10~ 5 for r? WMAP ), it is usually 
considered to be produced later through cosmic ray nu- 
cleosynthesis, but it is also known that the conventional 
processes are not sufficient for the observed abundance. 
This 6 Li problem was exacerbated by Q who found a rel- 
atively metallicity-independent abundance of 6 Li which 
is in contrast to the prediction of the cosmic ray synthe- 
sis scenario. In particular, a high abundance 6 Li/ 7 Li = 
0.046 ± 0.022 observed in the very metal-poor star LP 
815-43 with [Fe/H] = -2.74, is very hard to obtain. Un- 
conventional scenarios to enhance Li are investigated, 
for example, in [8( but they cannot solve the Li problem 
in any case. 

As has been discussed, the "Lithium problems", too 
much 7 Li and too little 6 Li produced in the standard 
scenario, do not have an astrophysical solution in a com- 
plete manner at present.* 1 In addition, the recent ob- 



# x Solutions by invoking not well measured reaction rates were 
investigated in fsL lldll . but without success at this stage. 
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servations seem to imply that stellar depletion of 7 Li is 
limited and that the metallicity dependence of 6 Li is only 
modest. These facts lead us to seek a solution to the 
Lithium problems by incorporating particle physics be- 
yond the standard model. Specifically, we reinvestigate 
the effects on BBN predictions of late-decaying parti- 
cles (LDPs) possessing a finite hadronic branching ratio. 
In the mid-eighties, Dimopolous, Esmailzadeh, Hall and 
Starkman (DEHS) showed that the mixed hadronic and 
electromagnetic decays of a massive particle at t ~ 10 5 s 
could reproduce the ratios of light elements EE EE] as 
then measured. In LDP nucleosynthesis, the products 
of particle decays scatter off the SBBN-produced light 
elements modifying the ALEs. Potentially, the modifi- 
cations in the ALEs could eliminate the existing incon- 
sistencies with current observations, since they predicted 
that a signature of such decays would be an anomalously 
high 6 Li to 7 Li ratio. 

LDPs, which we call here X, appear in widely con- 
sidered extensions to the standard model. For example, 
they are realised 

mis EE EE EE EE] m supergravity 

models where the next-to lightest sparticle (NLSP) de- 
cays into the lightest sparticle (LSP) with an extremely 
long lifetime (typically exceeding ~ 1 sec) , owing to 
the Planck-mass suppression of its interactions [l9l ]. In 
these theories the most favoured candidate for X would 
be a gravitino with a neutralino LSP [l6[, or a neu- 
tralino, a stau and a sneutrino as NLSPs with a grav- 
itino lsp [IE M, EE EE EE EE Ezt These are 

but a few of the plethora of subtly varying possibilities. 
Many have been discussed and their effects ri gor ously in- 
vestigated in studies conducted in the 1980 's Till [T2I . |28| . 
1990's [H, right up to recent times 0, ElT3!|7 with 
significant improvements being made at each stage (see 
also a text book [3l| and references therein, and a recent 
partial reconfirmation in Ref. [Hj]). Given the excellent 
agreement between SBBN and the measured (or inferred) 
abundances of D, 3 He and 4 He, the properties of LDPs 
are heavily constrained. Such constraints are valuable to 
theories beyond the standard model involving LDPs that 
are massive and weakly-interacting, which are difficult to 
study in collider experiments. 

The purpose of this paper is to utilise the state-of- 
the-art model-independent analysis of LDP nucleosyn- 
thesis performed by Kawasaki, Kohri and Moroi (KKM) 
[T3L [HI (incorporating important improvements of the 
original DEHS analysis), to identify the parameter space 
that is consistent with current observations. That is, we 
determine the range of values such as hadronic branch- 
ing ratio Bh, the lifetime tx, and the primordial energy 
density of X particles divided by primordial entropy den- 
sity Yx that solves the Lithium problems while leaving 
the abundances of H, D, 3 He and 4 He in agreement with 
observations. In addition we will reconfirm the original 
prediction of DEHS that a signature of the model is a 
high 6 Li/ Li ratio, and note that this prediction is now 
in agreement with observational data. Fundamental par- 
ticle physics models with LDPs would then need to lie 



within this allowed region of parameter space in order to 
solve the Lithium problems. Alternately, they might sat- 
isfy limits which have been presented elsewhere to ensure 
that they did not significantly change the ALEs, and look 
elsewhere for a solution of the Lithium problems. 

In Sec. HH we describe some key reactions in solving 
the Lithium problems. We carefully discuss how we can 
reduce uncertainties with regards the non-thermally pro- 
duced abundance of 7 Li and 7 Be resulting from the a-a 
collisions and make conservative predictions for the cor- 
responding effects on the 7 Li abundance. In Sec. IIII1 
we summarise recent observational data which we try to 
explain by BBN with LDPs. We present our results in 
Sec. IIVI and comment on the solution when some stellar 
depletion takes place in Sec. El We discuss our results in 
Sec. EH 



II. REACTIONS 

The reaction which is most significant in reducing the 
net 7 Li abundance is the neutron capture of 7 Be, 7 Be + 
n — > 7 Li + p, where the neutron is non-thermally pro- 
duced in the hadronic shower resulting from the hadronic 
decay of a LDP. Subsequently, 7 Li can be destroyed by 
thermal protons through the process 7 Li+p — > 2 4 He. For 
relatively high baryon to photon ratios, 77 > 3 x 10 -10 , 
the cosmic 7 Li abundance at much later times is mainly 
generated through electron capture by primordial 7 Be. 
This is because primordial 7 Li (but not 7 Be) is destroyed 
by the above thermal process involving p capture. How- 
ever, if large numbers of the non-thermal neutrons are 
emitted (e.g. by LDP decay) at the appropriate time, 
then the primordial 7 Be abundance can be converted to 
7 Li sufficiently early to be destroyed through thermal 
proton capture. This mechanism was originally identi- 
fied by Jedamzik in his pioneering work 15 1 for the first 
time and studied further in detail by [IE 1Z , EE HH (For 
other mechanisms to reduce the 7 Li abundance beyond 
the standard model, see also [H, HE HE HE HE El ) ■ 

The competing constraints on the properties of LDPs 
ultimately come from preventing the overproduction of 
D and 6 Li. These elements are also non-thermally pro- 
duced during the evolution of the LDP hadronic shower. 
D is directly produced through the spallation of 4 He by 
energetic nucleons, while 6 Li is mainly produced through 
the scattering of energetic (shower) T and 3 He off the 
background 4 He [ill El, EH- Such energetic T and 3 He 
are produced through the destruction of the background 
4 He within the LDP hadronic shower. 

Both 7 Li and 7 Be can also be non-thermally produced 
in a manner similar to 6 Li 

EH- 

If 4 He are produced dur- 
ing the evolution of the LDP hadronic shower through 
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He — » p/n + He +7rs (due to energetic shower 



nucleons), then energetic 4 He can collide with back- 
ground 4 He and produce 7 Li and 7 Be (and also 6 Li, al- 
though this is subdominant compared to the T - 4 He, and 
3 He- 4 He processes [HI EE, EH)- These non-thermal pro- 
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duction mechanisms were studied in detail in [13l . Il4j , in 
order to constrain the properties of the LDPs, and also 
applied to solve the 7 Li problem with some audacious ap- 
proximations in [l6j |. However, only if we adopt milder 
observational constraints on 7 Li (as in [H, c an wc 
ignore the intricacies of these production processes. 

Unfortunately, there is a severe lack of experimental 
data on the energy distribution of the 4 He in the final 
state of a-a inelastic scattering in the relevant energy 
regime. This energy distribution is essential in order to 
accurately compute the abundance of non-thermally pro- 
duced 7 Li and 7 Be. By making reasonable approxima- 
tions for the energy of the scattered 4 He, inferred from 
experimental data on similar processes involving colli- 
sions of energetic heavy ions, and utilising the theoretical 
properties of quantum chromodynamics (QCD), some of 
the current authors have investigated these non-thermal 
processes [3, [lfj- The constraints on the abundance 
of hadronically-decaying LDPs from predictions of the 
abundance of 7 Li and 7 Be were not significantly affected 
by these approximations, since the authors considered 
relatively generous observational constraints. 

However, for our purposes it is essential to more pre- 
cisely calculate the abundances of these elements. In par- 
ticular, a slight overestimation of the production of 7 Li 
and 7 Be may counterbalance their depletion by neutron 
capture, resulting in a worse fit to the more stringent data 
now available. Therefore we must reconsider the above 
ambiguities in order to better estimate the non-thermally 
produced 7 Li and 7 Be abundances. 

In the previous treatment [3, Hfl, there was a ten- 
dency for the estimated energy of the final state 4 He to be 
larger than the equipartition distribution in the center- 
of-mass (CM) system. This was because the authors 
used theoretical properties of quantum chromodynamics 
(QCD) to extrapolate results from high-energy experi- 
ments to lower energies until the extrapolation became 
kinematically inconsistent. However, whenever we in- 
clude non-standard processes, such as the aforementioned 
non-thermal a-a collisions, into the standard calcula- 
tions, it would be better to adopt most conservative ap- 
proximations that does not result in such an overestimate 
of the kinetic energy of the final state 4 He. Hence, we 
conservatively chose a (smaller) value of the 4 He kinetic 
energies between the experimentally-suggested QCD pre- 
diction [42j and the equipartition value in the CM sys- 
tem. For non-relativistic a-a collisions this tends to give 
smaller energies for the scattered 4 He, reducing the re- 
sulting abundance of Li and Be. To understand the sig- 
nificance of these processes, we also investigated scenarios 
without the a-a collisions. 



III. OBSERVATIONAL CONSTRAINTS ON 
LIGHT ELEMENT ABUNDANCES 

Our theoretical ALEs must be compared against the 
observational constraints on the abundances of D, 4 He, 



6 Li and 7 Li. The errors presented here are 1 a errors un- 
less otherwise stated. The subscripts "p" and "obs" refer 
to the primordial and observational values, respectively. 
The abundance of 3 He does not play any significant role 
in our conclusions. 

With regard to the mass fraction of 4 He, recently two 
groups reported new values of K, [IH Ej] by adopting 
quite new 4 He-emissivity data [45J. Izotov, Thuan and 
Stasinska reported two values, K[IZS1) = 0.2472±0.0012 
and F p (IZS2) = 0.2516±0.0011 [3 by using old and new 
He-emissivity data, respectively. Note that F p (IZS2) at 
face value excludes the SBBN prediction (~ 0.2484) even 
with various theoretical errors (~ 0.0004). We artificially 
incorporate a larger error into the value of F p (IZS2), 
to investigate conservative bounds, and call it "IZS3", 
Yp(IZS3) = 0.2516±0.0040, where the larger error 0.0040 
was adopted according to a discussion in [46[. The pos- 
sibility of these kinds of large errors was also discussed 
in gf . 

We adopt the following two deuterium abundances, 
Low (D/H) obs = (2.82 ± 0.26) x 10~ 5 as a most-recently 
reported value of the weighted mean [47| . and the 
more conservative value, High (D/H) b s = {3-98~to'^) x 

10" 5 m. 

As was discussed in Sec. HI we will compare 
both milder and more stringent constraints on 7 Li/H: 
logi ( 7 Li/H) obs = -9.63 ± 0.06 (MR) g|, and 
logi ( 7 Li/H) obs = -9.90 ± 0.09 (Bonifacio et al.) 

For the 6 Li abundance, we adopt 6 Li/ 7 Li = 0.046 ± 
0.022, which was recently observed in the very metal-poor 
star LP 815-43 with [Fe/H] = -2.74 0]. Again note that 
the value of 6 Li/ 7 Li calculated in SBBN (~ 3.3 x 10~ 5 ) 
does not agree with this constraint. 



IV. RESULTS 

Let us first summarise the basic framework of our 
study. We describe the properties of the LDP, which 
we call X, using only the following generalised parame- 
ters: -E v is, the (averaged) energy emitted in the form of 
visible particles (an invisible particle may also be emitted 
in some cases); Tx, the lifetime of X; B^, the branching 
ratio for X decay channels directly resulting in hadron 
production; the primordial abundance of X, which we pa- 
rameterise using the "yield variable" Yx = nx/s, which 
is defined at a cosmic time t <C tx- Here, nx is the num- 
ber density of X while s is the total entropy density. We 
know that the influence of unobservable decay products 
(e.g. neutrinos) has a negligible effect in our calculations 
(see [13, [H, [34j and references therein) . 

In our analysis, we calculate the primordial ALEs for a 
variety of different combinations of the above LDP model 
parameters, taking account of dissociation processes in- 
duced by the additional hadronic and electromagnetic 
interactions resulting from X decays, as discussed in 
[HI EH 1 with some important modifications discussed in 
Sec. HU We assume that two jets are produced in each 
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on the abundances of 4 He(ITS3), D, 7 Li and 6 Li, where 
in (a) we omit the a-a collisions as discussed in Sec. [Ill 
whereas in (b) we include the a-a collisions. Two com- 
binations of D and 7 Li constraints are displayed: one 
using Low (D/H) and 7 Li/H (MR), and a second using 
High (D/H) and 7 Li/H (Bonifacio et al.). Note that only 
upper bounds on B^E^Yx are provided by the K p (ITS3) 
and D/H contours. We also note that, even if we adopted 
other observational bounds on 4 He, such as ITS1 and the 
values published in [43|, [4(| , our results on the consistent 
region of parameter space would not be affected. 

In both FigfT] (a) and (b) (i.e. with and without 
the incorporation of a-a scattering), we identify a re- 
gion of parameter space that agrees with all obser- 
vational constraints, including those relating to 7 Li/H 
and 6 Li/H. Even if we adopted the more stringent 
constraints claimed by Bonifacio et al, there still re- 
mains such an allowed region of parameter space if we 
adopt the High(D/H) result. Any differences in these 
observationally-consistent regions of parameter space in 
figures (a) and (b) should be interpreted as due to theo- 
retical uncertainties in the hadron shower physics. For- 
tunately, such differences are minute. 

The 6 Li/ 7 Li contour located at smaller tx corresponds 
to the 2a lower bound 6 Li/ 7 Li > 0.002. It has the un- 
usual feature of being nearly perpendicular to the tx- 
axis. This is because the non-thermally produced 6 Li 
is uniformly destroyed by the standard thermal process 
6 Li(p, 3 He) 4 He, whose reaction rate is a steeply falling 
function of cosmic temperature. 

If we adopted a less-stringent lower bound, say 6 Li/ 7 Li 
> 6 x 10~ 5 (note that SBBN predicts 6 Li/ 7 Li ~ 3.3 x 10" 5 
for i] = 6.1 x IO" 10 ), it can explain Yp(ITS2). In turn we 
may say that we require models where 6 Li/ 7 Li > 6 x 10 -5 
in order to resolve the discrepancy between the Yp(ITS2) 
measurement and the corresponding SBBN prediction. 



FIG. 1: Bounds on BhE v i s Yx as a function of tx- In the 
panel (a), a-a collisions are excluded and in (b), a-a colli- 
sions are included. In each panel, the lower (upper) shaded 
region indicates the parameter space which is consistent with 
all observations when we adopt Low (D/H) and 7 Li/H (MR) 
(High (D/H) and 7 Li/H (Bonifacio et al.) ). The element's 
name is written beside each contour. 

hadronic decay of X, each with an energy £j e t = mj/2, 
with the hadronic branching ratio of X equal to B^. Here 
we set the mass of X to mx = 1 TcV and the energy 
converted into visible particles to E v i s — 2£j C t = mx- 
However, we conveniently find that even if we change 
mx, -Evis or £>h, the constraints on Bh-Evis^x are not 
significantly altered for tx ^ 10 6 sec. By comparing 
the results with the observations, we derived stringent 
constraints on both BhE v - ls Yx and tx- 

The shaded region indicates the part of parameter 
space that is consistent with observational 2a constraints 



V. STELLAR DEPLETION 

There is a possibility of late-time destruction of 6 Li 
and 7 Li in stars. Although we have argued in Sec. U 
that the detections of 6 Li in [7] show there is no stellar 
depletion, it can change for BBN with LDPs since 6 Li can 
be produced so much relative to 7 Li that the observed 
ratio 6 Li/ 7 Li is reproduced even though the depletion 
takes place. In fact, we find an interesting solution for 
which the constraints on the abundance of Li claimed 
by Bonifacio et al. and Ryan et al. are compatible with 
the Deuterium abundance constraint, Low (D/H), with 
some stellar depletion. 

According to 0, H§], rotational mixing within halo 
stars results in the depletion of Li, which can be paramc- 
terised by a depletion factor D7. The currently-observed 
7 Li abundance is then related to its primordial value by 
log 10 ( 7 Li/H) p = log 10 ( 7 Li/H) obs + D 7 . Because 6 Li is 
more easily dissociated than 7 Li, the depletion factor as- 
sociated with 6 Li, As ~ 2.5L> 7 [|. Thus ( 6 Li/ 7 Li) p = 
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FIG. 2: Results assuming stellar depletion. Here we have 
adopted the depletion factor Dj — 0, 0.1, and 0.2 (thin line, 
moderate line, and thick line). We see that the most severe 
bounds on 7 Li agrees with the theory with D? — 0.2. The 
star symbol means the allowed region. Here we have included 
a + a — » Li and Be processes. 

( 6 Li/ 7 Li) obs x lO 1 - 5 ^. 

In Fig. [21 we plot results for the ALEs corresponding 
to Dj = 0, 0.1, and 0.2. We see that the most severe 
experimental bounds on the abundance of 7 Li/H from 
Bonifacio et al. agree with the theoretical model for Dj > 
0.2, even if we adopt the Low (D/H) constraint. 

VI. DISCUSSION 

We have investigated a scenario where late-decaying 
particles (LDPs) of lifetime tx ^ 1 sec, possessing a fi- 
nite branching ratio into hadrons, B^, and emitting vis- 
ible particles of energy -E v is, are incorporated into the 
standard cosmological model. The scenario reaffirms 
the earlier prediction of DEHS that, as now observed, 
the 6 Li/ Li abundance ratio is much larger than the 
SBBN predicted value. We find that this allows us to 
solve the so-called "Lithium problems" , while simultane- 
ously remaining consistent with all other observational 
constraints on the abundances of light elements (ALEs) 
(within their 2 a error ranges). In other words, there 
exists a region of parameter space in which even the 



most severe observational bounds can be satisfied. It 
lies within the range 1.5 x 10~ 14 < B h (E vis /l GeV)Yx < 
3.0 x 10 -13 (where again Yx is the ratio of the primor- 
dial number density of LDPs divided by the entropy den- 
sity) and 3.0 < log 10 (tx /l sec) < 3.8 (with narrower 
ranges corresponding to observational constraints that 
are stronger than the most mild bounds discussed in this 
study). We derived the allowed regions in terms of the 
generalised parameters describing LDPs and converting 
them to some specific particle physics model parameters 
should be straightforward. This also means that the cur- 
rent study has independently confirmed the pioneering 
works flla . QjJ by Jedamzik and his collaborators, with 
an improved treatment of some of the key reactions in- 
volved, as discussed in §,H. 

The scenario has the possibility of the nonthermal- 
production of cold/warm dark matter as one of the decay 
products of the X, which must be important for forma- 
tions of small scale structures [l7|, [50| . 

Finally, since we have a lower limit on B^E^Yx to 
solve the Lithium problem, together with the constraint 
that the additional component of the energy density 
of the parent massive particle X not to produce too 
much 4 He, we can derive a lower limit on the hadronic 
branching ratio as B^ > 10~ 8 . The precise constraint 
Bh > 10~ 8 originates from the relatively reasonable as- 
sumption that we allow one additional neutrino species 
which contributes (approximately 15% of the total) to the 
energy density at the freezeout temperature of nucleons 
T = 0.1 MeV . 
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